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Axions and saxions from the primordial supersymmetric plasma 
and extra radiation signatures 
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We calculate the rate for thermal production of axions and saxions via scattering of quarks, gluons, 
squarks, and gluinos in the primordial supersymmetric plasma. Systematic field theoretical methods 
such as hard thermal loop resummation are applied to obtain a finite result in a gauge- invariant 
way that is consistent to leading order in the strong gauge coupling. We calculate the thermally 
produced yield and the decoupling temperature for both axions and saxions. For the generic case 
in which saxion decays into axions are possible, the emitted axions can constitute extra radiation 
already prior to big bang nucleosynthesis and well thereafter. We update associated limits imposed 
by recent studies of the primordial helium-4 abundance and by precision cosmology of the cosmic 
microwave background and large scale structure. We show that the trend towards extra radiation 
seen in those studies can be explained by late decays of thermal saxions into axions and that 
upcoming Planck results will probe supersymmetric axion models with unprecedented sensitivity. 

PACS numbers: 14.80. Va, 11.30.Pb, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

There are several hints towards physics beyond the 
standard model (SM). One of them is the strong CP prob- 
lem. If this problem is solved via the Peccei-Quinn (PQ) 
mechanism, the axion a arises as the pseudo-Nambu- 
Goldstone boson associated with the U(l)pQSymmetry 
broken spontaneously at the PQ scale /pq [J • An- 
other attractive extension of the SM is supersymmetry 
(SUSY) 0-01. In conceivable settings with both the PQ 
mechanism and SUSY, the pseudo-scalar axion is part 
of a supermultiplet in which also its scalar partner, the 
saxion a, and its fermionic partner, the axino a, appear. 
The energy density of the early Universe can then re- 
ceive contributions from coherent oscillations of the axion 
field [H, [H, 0] and the saxion field [§-[l2j and from thermal 
production of axions [l], [l3|-[l7j , saxions @, [l(| EH , and 
axinos [l9l - [29| in the hot primordial plasma. 

Here we calculate for the first time the thermal produc- 
tion rate of axions and saxions via scattering processes of 
quarks, gluons, squarks, and gluinos in a gauge-invariant 
way consistent to leading order in the strong coupling 
constant g s . In our calculation we use hard thermal loop 
(HTL) resummation 30] and the Braaten-Yuan prescrip- 
tion [3l| to account systematically for screening effects 
in the quark-gluon-squark-gluino plasma (QGSGP). This 
method was introduced on the example of axion pro- 
duction in a hot QED plasma (3l|; see also Ref. [HJ. 
Moreover, it has been applied to calculate the thermal 
production of gravitinos }33 - l35j and axinos [HJ in SUSY 
settings and of axions in a non-SUSY quark-gluon plasma 

(QGP) Ql M- 

Based on our result for the thermal axion/saxion pro- 
duction rate, we determine the respective thermally pro- 
duced yields and estimate the decoupling temperature of 
axions and saxions from the thermal bath. While both 
axions and axinos are promising dark matter candidates 
(cf. [H, H, H3, Hl| and references therein), saxions can be 



late decaying particles with potentially severe cosmolog- 
ical implications. For example, energetic hadrons and 
photons from saxion decays during or after big bang nu- 
cleosynthesis (BBN) can change the abundances of the 
primordial light elements [i"lj ]. Moreover, photons from 
saxion decays can affect the black body spectrum of the 
cosmic microwave background (CMB) for a saxion life- 
time of 10 6 s < r a < 10 13 s [ill HI or may contribute 
either to the diffuse Y(7)-ray background or as an addi- 
tional source of reionization for r a > 10 13 s (ill. I39l44li. 



In scenarios in which the decay mode into axions is not 
the dominant one, saxion decays ma y also p roduce signif- 
icant amounts of entropy 0, [H, 1391.143-144 1 . This can di- 
lute relic densities of species decoupled from the plasma 
and also the baryon asymmetry r/. Then, r a < 1 s is 
imposed by successful BBN which requires a standard 
thermal history for temperatures below T ~ 1 MeV. 

In this work, however, we look at scenarios in which 
saxions (from thermal processes) decay predominantly 
into axions. Moreover, we still focus on decays prior 
to BBN and compute the additional radiation provided 
in the form of the emitted relativistic axions. Such a 
non-standard contribution AN c g to the effective num- 
ber of light neutrino species N c g from decays of ther- 
mal saxions into axions was previously considered in 
Refs. [1, [ll|, HH Ell- Applying our new result for the 
thermally produced saxion yield and new cosmological 
constraints on AJV c ff imposed by recent studies of BBN, 
the CMB, and large scale structure (LSS) [iH4j|. we 
present updated limits on the PQ scale /pq, the saxion 
mass m a , and the reheating temperature Tr after infla- 
tion. 

Interestingly, precision cosmology (48145211 and recent 
studies of the primordial 4 He abundance |46l . [47j show a 
trend towards a radiation content that exceeds the pre- 
dictions of the SM. In fact, such an excess can be ex- 
plained by the considered saxion decays into axions. The 
observed trend may thus be a hint for the existence of a 
SUSY axion model. Here results from the Planck satel- 
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lite mission will be extremely valuable, which will come 
with an unprecedented sensitivity to the amount of extra 
radiation at times much later than those at which BBN 
probes this quantity. Based on a forecasted 68% confi- 
dence level (CL) sensitivity of AN e g = 0.26 [53j, [5J], we 
indicate parameter regions of SUSY axion models that 
will be tested by results from the Planck satellite mis- 
sion expected to be published in the near future. 

The remainder of this paper is organized as follows. 
In the next section we consider interactions of the PQ 
supermultiplet and decay widths for saxion decays. In 
Sects. Hill and II Vl our calculations of the thermal produc- 
tion rates of saxions and axions are presented. We com- 
pute the associated yields in Sect. IVl and use the results 
to estimate the saxion/axion decoupling temperature in 
Sect. IVI1 Then we explore AN e s provided in the form of 
axions from saxion decays and possible manifestations in 
studies of BBN and of the CMB and LSS. Here we com- 
ment on potential restrictions which can emerge from 
overly efficient thermal gravitino/axino production and 
describe exemplary settings that allow for a high reheat- 
ing temperature of T R ~ 10 8 -10 10 GeV. In Sect. NTU\ 
we compare the relic density of axions from the misalign- 
ment mechanism with the ones of thermal axions and of 
non-thermal axions from saxion decays. Our conclusions 
are given in Sect. |IX] 



II. PARTICLE PHYSICS SETTING 

In a SUSY framework, the U(1)pq symmetry is ex- 
tended to a symmetry of the (holomorphic) superpoten- 



component spinor notation: 



tial and thereby to its complex form U(l)pq (55j. In 
the case of unbroken SUSY, this implies the existence 
of a flat direction and thereby a massless saxion field. 
Once SUSY is broken, this flat direction gets lifted, re- 
sulting in a model-dependent mass of the saxion m CT . 
For example, m is expected to be of the order of the 
gravitino mass fn?.ip. in gravity-mediated SUSY breaking 
models [HI [H, Here we do not look at a specific 
model but treat m a as a free parameter. 

In this work we consider the particle content of the 
minimal supersymmetric SM (MSSM) extended by the 
PQ superfield A = (a + m)/V2 + y/26a + F A 69, where 
8 denotes the corresponding fermionic superspace co- 
ordinate and Fa the chiral auxiliary field. The in- 
teractions of A with the color-field-strength superfield 



W 



~9 b - 



D b 9 - a^OG" + mo^DyS* are given by the 



effective Lagrangian 



V2c 



PQ 



8tt/pq 



d 2 9A W b W b + h.c. 



(1) 



/•int _ a s 

LpQ ~8^ PQ 



7.h 



a {G b ^G% - 2D»D» - 2^ 7 ^mJ 



a(G b ^G b llv +2f M ^D ll ~g b M ) 



2a M D b g b M 



(2) 



where G b 
with a 



.,„. e» Vfm G b >>°/2 and D b = V.^/J VS/, 
sum over all squark fields q and the SU(3) C 
generators T b j in their fundamental representation; the 
subscript M indicates 4-component Majorana spinors. 
Note that we use the space-time metric g M „ = g^ v = 
diag(+l, — 1, — 1, — 1) and other conventions and nota- 
tions of Ref. [6j and - except for a different sign of the 



0123 



T - of Ref. [3] • For compar- 



Levi-Civita tensor e 1 
isons with similar £pq expressions given in Refs. [25l.l28j. 
we remark that the second term in the brackets in the sec- 
ond line of ([I]) can be written as Dy(g b M j ,1 j 5 g b M ). How- 
ever, our result for the saxion-gluino- interaction term dif- 
fers from the corresponding terms in [25j and [28j by fac- 
tors of —2 and —1, respectively. Moreover, our findings 
for the axino interactions differ by a factor of — I from the 
ones in [25l . |28| . This may result partially from metric 
conventions: If we translate ^ into the corresponding 
expression valid for g^ v = g^ = diag(— 1, +1, +1, +1) 
using Appendix A of Ref. [6j , the sign of our result for the 
axino-gluino-gluon-interaction term will change, whereas 
all other terms in ([2]) will not be affected. 

In the following we focus on hadronic or KSVZ axion 
models [5(1 [57j in a SUSY setting in which the effec- 
tive Lagrangian describes the relevant interactions 
even in a conceivable very hot early stage of the pri- 
mordial plasma with temperatures T not too far below 
/pq. Note that we do not consider scenarios with a 
radiation-dominated epoch with T above the masses of 
the heavy KSVZ (s)quarks tUq q such as those consid- 
ered in Ref. 

Next we address interactions between axions and sax- 
ions in models with N SM-gauge singlet PQ multiplets 
$i with PQ charges qi and vacuum expectation values 
(VEVs) Vi after PQ symmetry breaking. Near the VEVs 
(after integrating out the heavy fields in the PQ sector), 
the scalar parts of $i can be parametrized as 



Vi exp 



qi(a + ia) 



V2' 



VPQ 



(3) 



where b is a color index, g the gluino field, D the 
real color-gauge auxiliary field, the gluon-field- 

strength tensor, D„ the corresponding color-gauge co- 
variant derivative, and a s = gg/(47r). After performing 
the integration, we get for the propagating fields in 4- 



Here the canonical PQ charge normalization requires 
qf = 1 for the smallest qi and wpq = v 'i Qi results 

from the requirement of canonically normalized kinetic 
terms for the axion and the saxion. Moreover, from the 
kinetic terms of the PQ fields, one finds that interactions 
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between axions and saxions can emerge as follows [2 11 ] 



JY 



\[2x ^ 



i=l 



VPQ 



(4) 



where x = JV qf v~ /vpQ. The strength of these in- 
teractions thus depends on the model. For example, 
x — (v\ — v%)/v P q in models whose superpotentials con- 
tain the term /ci?($i3>2 — v P q/2) with a Yukawa coupling 
k, two PQ fields with qi = —q% = 1 and a SM-gauge sin- 
glet field R with q R = 0. This illustrates that x <C 1 is 
possible if v x ~ v 2 ~ wpq/v^ [U HI HI]. On the other 
hand, in a KSVZ axion model with just one PQ scalar 
(with v = vpq and q — 1) (Toj . one finds x = 1, which is 
the value that we will consider in Sect. IVIII below. 

Let us now relate the scale Upq, imposed by canon- 
ically normalized kinetic terms, to /pq, defined by the 
form of the prefactor of the effective axion-gluon interac- 
tions in In a KSVZ model, those interactions emerge 
from axion couplings to the heavy KSVZ quarks which 
are described by contributions to the superpotential of 
the form h^iQ^Qn with a Yukawa coupling h, q\ = 1 
and heavy quark multiplets Ql and Qr with color charge 
and PQ charges qQ = —1/2. Considering the resulting 
Lagrangian that describes the interactions of </>i with the 
KSVZ quarks, one sees that the spontaneous breaking of 
the PQ symmetry results in heavy Dirac KSVZ quarks 
with a mass uiq = hv\. Integrating out loops of such 
heavy quarks, one finds the effective Lagrangian describ- 
ing axion interactions with gluons 



C 



int 



hv\a s 



8irV2 



mQVpQ 



(5) 



where the KSVZ quarks have been assumed to be in 
the fundamental representation of SU(3) C . For /pq = 
V^vpq, one thus recovers the well-known form of the 
corresponding interaction term as given in 

Numerous laborat ory , astrophysical, and cosmological 
studies point to [53, [6(| 



/pq > 6 x 10 8 GeV. 



(6) 



This corresponds to an upper limit of about 10 meV on 
the axion mass, 



m a — 6 meV 



10 9 GeV \ 
/pq / 



(7) 



1 Here we focus on N, 



Q 



and Q R . For Nq > 1. 



1 heavy KSVZ (s)quark multiplets Ql 
/pq -> /PQ/Wq, e.g., in 0, ©, ©, 



and |(2J in line with an additional factor of Nq on the right-hand 
side of |(5}. Using this /pq definition, there are no modifications 
of the relation /pq = \72"pq and of JSJl below for Nq > 1. 



and implies that axions are stable on cosmological 
timescales. Because of the larger mass of the saxion, 
its lifetime r CT is typically smaller than the age of the 
Universe and governed by the following decay widths. 
From ffl one obtains the width for the saxion decay into 



axions J 



64t™ pq 



327T/2 



(8) 



PQ 



and from ([2]) the width for the saxion decay into gluons, 



>39 



167T3/IV 



(9) 



For KSVZ fields that carry an non-zero electrical charge 
eQe with e = y/Ana and the fine-structure constant a, 
the saxion can decay into photons via KSVZ quark loops. 
After integrating out those loops, we find the associated 
width 



a— ^77 



9eQa 2 ml 

647T 3 /, 2 



(10) 



PQ 



If x > 0.2, the saxion decay into axions governs T ai 
which is the case on which we focus in this work. Indeed, 
in the region m a > 1 GeV in which the competing decay 
a — > gg is possible, such x values imply the branching 
ratio BR((T — > aa) > 0.9. For m a below the threshold to 
form hadrons, where a — > 77 is the competing decay, the 
decay into axions governs r CT for even smaller values of x, 
e.g., for eg = 1 and x = 0.02, we still find the branching 
ratio BR(er -> aa) > 0.9. 



III. THERMAL SAXION PRODUCTION 

Let us now calculate the thermal production of saxions 
in the primordial SUSY QCD plasma. Assuming that in- 
flation has governed the earliest moments of the Universe, 
any initial population of saxions has been diluted away 
by the exponential expansion during the slow-roll phase. 
After completion of the reheating phase that leads to a 
radiation-dominated epoch with an initial temperature 
Tr, the thermal production of saxions starts to become 
efficient. In fact, we focus on cosmological settings in 
which radiation governs the energy density of the Uni- 
verse as long as this production mechanism is efficient 
(i.e., for T down to at least T ~ 0.01 Tr). While inflation 
models can point to Tr well above 10 10 GeV, we consider 
the case Tp < /pq such that no PQ symmetry restora- 
tion takes place after inflation. Moreover, Tp < tuq q 
is assumed in line with our comments on the considered 
KSVZ axion model settings in the previous section. 



2 Our result ((Sj) agrees with the ones of Refs. fiol. 1 3S|| . where x ■■ 
and /pq = tipq, and of Refs. [HI. [6lll63l . where F a = fpQ- 
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The calculation of the thermal production of saxions 
with E > T follows closely [l7|, where thermal axion 
production in a SM QGP is considered. From ((2J we get 
the relevant 2 — > 2 processes shown in Fig. [TJ Additional 
processes exist but can be accounted for by multiplying 
the squared matrix elements of the shown processes with 
appropriate multiplicity factors. The squared matrix ele- 
ments of the shown processes are listed in Table HI where 
s = {Pi +P 2 f and t = (Pi-Pzf with Pi, P 2 , P 3 , and P 
referring to the particles in the given order. Working in 
the limit, T>m,, the masses rrii of all MSSM particles 
involved have been neglected. Sums over initial and final 
spins have been performed. For quarks and squarks, the 
contribution of a single chirality is given. The obtained 
squared matrix elements can be calculated conveniently, 
e.g., with the help of FeynArts [63[ and FormCalc [Hj]. 

The results for processes A, C, E, and G given in Ta- 
bic U point to potential infrared (IR) divergences. Here 
scree ning effects of the plasma become relevant. In 
Ref. foCl loll] a systematic method is introduced to ac- 
count for such screening effects in a gauge-invariant way. 
Following Ref. [HI , we introduce a momentum scale fc cu t 
such that g s T k cut < T in the weak coupling limit 
g s <C 1. This separates soft gluons with momentum 
transfer of order g s T from hard gluons with momentum 
transfer of order T. By summing the respective soft and 
hard contributions, the finite rate for thermal production 
of saxions with E > T is obtained in leading order in g s , 



E- 



d 3 p 



E 



d 3 p 



E- 



soft 



d 3 p 



(11) 



hard 



which is independent of k cut - 

In the region where k < k cnt , we use the optical the- 




Process B: qi + q 3 ■ — > g a + a 




Process C: qi + g a — s> qj + a (Crossing of B) 
Process D: g a + g b ->■ g c + a 




g & O) g g <ep % s 

Process E: g a + g b — >• g c + a (Crossing of D) 
Process F: g 4 + q~j — > g" + a 



TABLE I. Squared matrix elements for saxion (er) production 
in 2-body processes involving MSSM quarks and squarks of 
a single chirality (qi, qi), gluons (g a ), and gluinos (g a ) in 
the high-temperature limit, T ^> rrii, with the SU(3) C color 
matrices f abc and Tj t . Sums over initial and final state spins 
have been performed. 



Process G: q~i + g a — >• qj + a (Crossing of F) 
Process H: q~i + q~j — > g a + a 



Label i 



Process i 



A 
B 
C 
D 
E 
F 
G 
H 
I 
.1 



9 +9 
qi + q 3 

qi + 9 a 

F+g" 

9 a +9 b 
qi + qj 

q>+g a 

qi + q~j 

g a + q> 
g a + q t 



g +° 

g a +a 
g c + a 



9 

g a 

9 a 

qj 
qj 



- (7 

- a 

- (j 

- (j 

- a 

- a 



A (s*+st+t~y 1 f abci2 
* st(s + t) \J I 

- +2t + s] IT"! 2 



2f -2s-t) \T^\ 2 
'^-+2t + s^ \f abc \ 2 
2d-2a-t) \f abc \ 2 



-2( t -+t)\T^ 



s\T^\ 2 
-t\T-\ 2 



Process I: g a + qi — > qj + a (Crossing of H) 
Process J: g a + q~i — !> qj + a (Crossing of H) 

FIG. 1. The 2 — > 2 processes for saxion production in a 
SUSY QCD plasma. Additional processes are included in 
terms of multiplicities in our calculation of the thermal pro- 
duction rate: Process C with antiquarks q~ij replacing qtj, 
process G with antisquarks q it j replacing q~ij, process H with 
antisquarks/quarks q~i/qj replacing q~i/qj, and processes I and 
J with qi and q~j replacing qi and q~j, respectively. 
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FIG. 2. The saxion self energy used to compute the leading 
contribution to the thermal production rate of hard saxions. 
The blob indicates the HTL-resummed gluon propagator. 



orem to obtain the soft contribution from the imaginary 
part of the saxion self energy shown in Fig. [5] Since only 
one gluon can carry a soft momentum, we need to use 
the HTL-resummed propagator only once. Using fc cut as 
the ultraviolet cutoff, we get 



E 



d 3 p 
= Ef B (E) 



(2tt) 3 



Imn ff (£ 

- 1)T 



,P)\ 



8192tt8/2 q 



m 9 



1.379 



(12) 
(13) 



where the squared SUSY thermal gluon mass is given by 
m 2 = g 2 T 2 (N c + n f )/6 for N c = 3 colors and n f = 6 light 
quark flavors and the equilibrium phase space density for 
bosons (fermions) by / B (f)(-E') = [exp(-E/T)q=l] -1 . More 
details on the way in which this calculation is performed 
can be found in Refs. [3l|, [32J, l35l. l36j. 

In the region where k > fc cu t we can use zero tem- 
perature Feynman rules since fc cu t provides an IR cut- 
off. From the matrix elements given in Table HI weighted 
with appropriate multiplicities, statistical factors, and 
phase space distributions, we get the (angle-averaged) 
hard contribution 



E 



d 3 p 



hard 



2(2tt) s 



47T 



d 3 P] 



n (2^2^ 



x (2tt) 4 S 4 {P 1 +B2-P3- P)&{k - fccut) 



Y J h{E 1 )f 2 {E 2 )[l ± f 3 (E 3 )]\M 1+2 ^ : 



3+ff| 



(14) 



E g 6 a (NZ-l){N c + n f ) I f B (E)T 3 



512tt 7 /2 
h{E)T 3 



V. 



32tt 

(i) , 

BBB ' 



PQ 
T 2 



48tt 



Ia( F -]+--27- 



'FBF 



v cut 

A3) 
'BBB 



17 

T 



f (3) ) 



C(2) 
»/(4bf + h 



(2) 

FFB^ 



Nr. 



n f 



with Euler's constant 7, Riemann's zeta function £(z), 



1 



BBB(FBF) " 32^3^2 J Q 



dE : . 



e+e 3 



dEx In 



\E\ — E 3 



E 3 



-Q(Ei - E 3 )—[f BBB{FBF) E 2 2 {El + El)] 



+ Q(Ea - Ei) -j-j^-[f BBB (F B F)E 2 (Ef + El)] 
+e(E-E 1 )^-[f BBB{FBF) (E 2 1 E 2 -E 2 E 2 )]\, (16) 



r(2) 



dE 3 



E+E 3 



dE 2 f F 



'FBF(FFB) - 96^3^2 J J JFBF(FFB) 

{6(E - E 3 ) [E 2 3 (E 3 - 3Ei) - G(E 2 - E){E 2 - E) 3 ] 
+ e{E 3 - E) [{E - 3E 2 ) E 2 + e(E - E 2 ){E 2 - E) 3 ] 
+ [Q(Ea - E 2 )G{E - E 3 ) - <d(E 3 - E)G{E 2 - E 3 )} 

x [(E 2 ~E 3 ) 2 (E 2 + 2E 3 )-3(E 2 2 -E 2 3 )E]}, (17) 



r(3) 



1 



BBB(FFB) 32^3^2 J q 



dE 3 



E+E 3 



dE 2 ft 



BBB(FFB) 



X *&{E-E 3 ) 



EM 



&(E 3 - E) 



E 2 E 2 



E 3 + E E 3 + E 

+ [&{E 3 - E)0{E 2 - E 3 ) - B{E - E 3 )B{E 3 - E 2 )} 

x {E 2 - E 3 ) [E 2 (E 3 -E)- E 3 (E 3 + £)]}, (18) 

/bbb,fbf,ffb = fi(E 1 )f 2 {E 2 )[l ± f 3 (E 3 )}. (19) 



The sum in (|14[) is over all saxion production processes 
1 + 2 — » 3 + a viable with ([2]) . The colored particles 1-3 
were in thermal equilibrium at the relevant times. Per- 
forming the calculation in the rest frame of the plasma, 
fi are thus described by f F i B depending on the respec- 
tive spins. Shorthand notation (fT9"|) indicates the corre- 
sponding combinations, where + (— ) accounts for Bose 
enhancement (Pauli blocking) when particle 3 is a boson 
(fermion). With any initial saxion population diluted 
away by inflation and for T well below the saxion de- 
coupling temperature Tb (which will be determined in 
Sect|V}, we can neglect saxion disappearance reactions 
and Bose enhancement by saxions, since the saxion phase 
space density f a <C f F / B and 1 + f a s» 1. 



IV. THERMAL AXION PRODUCTION 

The calculation of thermal axion production in the pri- 
mordial SUSY QCD plasma proceeds analogously to the 
saxion calculation presented in the previous section. Af- 
ter substituting the saxion a by the axion a, the Feynman 
diagrams can be read directly from Figs. Q] and [2] with one 
modification: there is no gluino-gluino-gluon-axion ver- 
tex and thus no quartic interaction such as the one that 
contributes to processes D and E in the saxion case. 

Although the Feynman rules for the axion interactions 
derived from (j5J| differ from the ones describing saxion in- 



teractions, we obtain squared matrix elements for the ax- 
ion production processes in the high-temperature limit, 
T rrii, that agree with the ones for the correspond- 
ing saxion production processes given in Table HI More- 
over, we find that both the soft and the hard contribu- 
tions to the thermal production rate of hard axions agree 
with (fT3| and (|T5|) . respectively. Our result for the ther- 
mal axion production rate E dW a /d 3 p thus agrees with 
the one for the thermal saxion production rate obtained 
above. This implies an agreement of the associated ther- 
mally produced yields of axions and saxions prior to de- 
cay, which will be calculated in the next section. 

Before proceeding let us stress that we can neglect pro- 
duction processes like tttt — > na in the primordial hot 
hadronic gas [3, EH because of the /pq limit ((6]) . More- 
over, Primakoff processes such as e~7 — > e~a are not 
taken into account since they are usually far less efficient 
in the early Universe [l3| . 



V. THERMAL SAXION/ AXION YIELD 

Let us now calculate the thermally produced (TP) sax- 
ion yield Yj p = n a /s, where n a is the corresponding sax- 
ion number density and s the entropy density. With the 
results obtained in the two previous sections, we know 
beforehand that this yield prior to decay agrees with the 
thermally produced axion yield Y" a TP = n a /s. While the 
calculation and results are indeed valid for both saxion 
and axion, we focus on the saxion case. 

For T sufficiently below the saxion decoupling temper- 
ature Td, the evolution of the thermally produced n a 
with cosmic time t is governed by the Boltzmann equa- 
tion 
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dn a 
~dt 



+ ZHn 



d 3 p 



dW a 
d 3 p 



(20) 



Here H is the Hubble expansion rate, and the collision 
term is the integrated thermal production rate: 



W„ = 



9C(3)g, 6 T 6 

256tt 7 /2 



PQ 



hi | ^— ) + 0.4:505 

m 9 



(21) 



Assuming conservation of entropy per comoving volume 
element, (|20|) can be written as dYj p jdt — W a /s. Since 
thermal saxion production is efficient only in the hot ra- 
diation dominated epoch with temperatures well above 
the one of radiation- matter equality, T mat=ra d, we can 
change variables from cosmic time t to temperature T 
accordingly. With an initial temperature Tr at which 
YJ p (Tr) ~ 0, the relic saxion yield prior to decay is 



Y. 



TP 



1.33 x 10" 3 5 f In 



To, 

1.01 

9s 



dT 



W a {T) 
Ts{T)H(T) 

10 n GeV \ 2 

/pq / 



T, 



R 



10 8 GeV 




(22) 



10" 

Tr [GeV] 

FIG. 3. The relic saxion yield prior to decay originating from 
thermal processes in the primordial plasma for cosmological 
scenarios characterized by different Tr values covering the 
range from 10 6 to 10 10 GeV. The dash-dotted, dashed, and 
solid lines are obtained for / PQ = 10 10 , 10 11 , and 10 12 GeV. 
The relic axion yield Ya q ^ TP from thermal processes agrees 
with Y^ TP and can thus be read from this figure as well. 



with a fiducial temperature T\ ow well below Tr and well 
above T a , which we use to denote the temperature of 
the primordial plasma at t — r CT : T a <C T ow <C Tr. In 
the case of the axion, Tj ow = T mat=ra( j can be used since 
its lifetime exceeds the time of radiation-matter equality 
significantly. Note that the resulting saxion/axion yield 
is insensitive to the exact choice of T\ ow for T\ ow < 0.01 Tr 
since additional contributions from thermal production 
at T < 0.01 Tr are found to be negligible. 

Figure [3] shows the saxion yield (|22| for /pq = 10 10 , 
10 11 , and 10 12 GeV as the diagonal dash-dotted, dashed, 
and solid lines, respectively. Here we compute (|2"2")l 
with g s = <7 s (Tr) = y/ 47tq; s (Tr) evaluated according 
to its 1-loop renormalization group running within the 
MSSM from a s (mz) — 0.1176 at the Z-boson mass 
mi = 91.1876 GeV. The applied methods [3(||3H require 
g s < 1, so that ([22]) is most reliable for T R > 10 6 GeV. 
For lower Tr values such that g s (Tn) > 1, one encoun- 
ters an artificial suppression of Yj p and even unphysi- 
cal negative values, which can be seen directly from the 
logarithmic factor in (|22l) . This is a well-known limita- 
tion of this technique (cf. [13, HU) that calls for gener- 
alizations of the gauge-invariant methods introduced in 
Refs. [13, HH modified to extend the applicability beyond 
the weak coupling limit. 

Note that ([22]) is only valid if T R < T D , because oth- 
erwise saxion annihilation processes neglected in (|14|) are 
important. For Tr > Tp> saxions were in thermal equilib- 
rium in the early Universe before decoupling as thermal 
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relics. In fact, for m a <C 7b, saxions decouple as a rela- 
tivistic species. The yield is then given by 



yeq _ 



1.2 x 10 



-3 



(23) 



as indicated by the horizontal lines in Fig. [3J The yield 
from thermal production cannot exceed the equilibrium 
yield, so that (1231) represents an upper limit. In scenar- 
ios with values of Tr for which ([22]) is equal or larger 
than (|23p . saxion disappearance reactions have to be 
taken into account. The resulting yield would show 
almost the same Tr dependence as the one in Fig. [31 
but with a smooth transition in the Tr range in which 

yTP ^ T^oq 

For axions, m a <C 7b is always satisfied so that they 
will be hot thermal relics with Y° q = Y° q as given 

in d23[) if T R > T D . In fact, the relic axion yield Ya° q/TP 
from thermal processes agrees fully with FCT q ^ TP described 
above and can be read from Fig. [3] as well. 



VI. DECOUPLING TEMPERATURE 

Considering Fig. [31 one finds that the kinks indicate 
critical Tr values. For a given /pq, the associated critical 
Tr value separates scenarios with thermal relic saxions 
from those in which saxions have never been in thermal 
equilibrium. We thus use the positions of the kinks as 
/pQ-dependent estimates of the saxion decoupling tem- 
perature. Our numerical results are well described by 



7b « 1.4 x 10 9 GeV 



fpQ 



10 11 GeV 



(24) 



This is similar to the estimate of the axino decoupling 
temperature in Ref. [l9j . Such an agreement was ex- 
pected and used to provide estimates of the thermally 

" HI. An- 



produced saxion yield Y a in Rcfs. [lfl [11|, [22|, 
other recent study applies the thermally produced axino 
yield Y a TP Q to estimate Yj p ~ (2/3) F a TP [H3|0 With 
our results illustrated in Fig.[3]above, one can now see ex- 
plicitly the similarity between Yj p and the correspond- 
ing axino yield F a TP illustrated in Fig. 4 of Ref. [2J] . 

In light of Sect. IIV1 it is clear that ([2~4"[) describes the 
axion decoupling temperature in the considered SUSY 
settings as well. When comparing ([24]) with the axion 
decoupling temperature in non-SUSY scenarios, given in 
(15) of Ref. [17|, we find only small differences. In 
1, the additional diagrams in the 



Eq, 

fact, for a fixed g 
SUSY case (which lead to a different thermal gluon mass 



3 Note that / P Q in [Io ll39| and F a in [ll|, l62ll correspond to our 
"PQ = /pq/v2 and thereby differ by l/y/2 from our /pq. With 
these differences in the definitions of the PQ scale, we find that 
the Yj p estimates in Refs. [TTj, EJ [H EI exceed the result (22} 
of our calculation by about a factor of two for fixed /pq and Tr. 



m g also) increase the collision term for thermal axion 
production W a only by at most 30% with respect to 
Eq. (12) of Ref. [H obtained for the non-SUSY case. 
Note also that both Yj p and Y a eq are normalized to an 
entropy density s(T) which is more than two times larger 
in the SUSY case than in the non-SUSY case due to the 
additional sparticles which can all be considered to be 
relativistic at very high temperatures such as the axion 
decoupling temperature. 



VII. ADDITIONAL RADIATION FROM 
SAXION DECAYS 

As already mentioned in the Introduction, axions from 
late saxion decays can provide additional radiation al- 
ready prior to BBN and later on as well. The amount 
of additional radiation is usually expressed in terms of 
a non-standard contribution AN e g to the effective num- 
ber of light thermally excited neutrino species N e g. It is 
defined in relation to the total relativistic energy density 



Prad(T) 



1 



p.,(T) (25) 



with the photon energy density p 7 and the temper- 
atures of neutrinos T v and of photons T. At T > 
1 MeV (before neutrino decoupling and e + e~ annihila- 
tion), N cS = 3 + AiVeff and T v = T. These relations 
change to T u = (11/4)~ 1//3 T after neutrino decoupling 
and to iVcff = 3.046 + AiVeff [65[ because of residual neu- 
trino heating by e + e~ annihilation. 

At a photon temperature T <T a , the energy density of 
relativistic non-thermally produced (NTP) axions from 
saxion decays Pa TP (T) yields 



AN cS (T) 



120 



NTP 
Pa 



(T). 



(26) 



Working in the sudden decay approximation, all ther- 
mally produced saxions are considered to decay instan- 
taneously at t T(j (where T = T a ). If the saxions are 
non-relativistic when decaying dominantly into two ax- 
ions, the initial axion momentum is p a (T a ) — m a /2 and 

1/3 rp 

(27) 



P NTP (T) = ^ 



9*s(T) 



g*s(T a ) 



' 9*s{T) ' 


4/3 


r) 


_g*s(T a )_ 


( 


\T a J 



— n eq / TP 



p^ ip (T ff ) (28) 



with p^ /TP (T CT ) = m a Y^ /TP s(T a ) and Y a NTP = 2Y ff cq/TP . 
Here g*s denotes the number of effectively massless de- 
grees of freedom such that s — 27r 2 (? >t sT 3 /45. For r CT ~ 
1 /T a ^ aa given by ([8]) and with the time-temperature re- 
lation in the radiation-dominated epoch, we obtain 



T a ~ 10.6 



MeV( 



Vl GeV 



/2/ 10 10 GeV \ 
V fpQ/* J 



10.75 



-1 1/4 



9*{Ta) 



(29) 
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and 



AA cff (T) 



0.95 ^1 00 GeV \ 1/2 
m a ) 




L 
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'g*s{T)~ 


4/3 9*{t^ 




10.75 _ 





(30) 



where is the effective number of relativistic degrees of 
freedom governing the energy density. 

Focussing on saxions from thermal processes, the max- 
imum AiVeff emerges for scenarios with Tr above the 
decoupling temperature so that the thermal relic 
yield (j2"3")l applies: 



AA off (T) 



1.14 ^100 GeV\ 1/2 
m a ) 



PQ 



10 11 GeV 



(f 
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g*s{T)~ 
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1/4 



(31) 
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For T R < T D on the other hand, the yield ((221) leads to: 



AiVeff(T) 



12.6 ffs 6 ln(Mi 



1 GeV\ 1/2 /10 n GeV\ 
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(32) 



Thermal relic saxions are non-relativistic when decay- 
ing if their average momentum at T a satisfies 



(p(T ff )> = (p(T D )) 



g*s{Tg) 
g*s(Tu) 



1/3 



(33) 



Since those saxions decouple as a relativistic species (pro- 
vided m a <C T D ) at a very high temperature (|2~4"|) with 
a thermal spectrum, (p(Tb)) = 2.701 Td and <?*s(Tb) — 
232.50 Using (|29l) . we can express (|33l) in terms of the 
following m^-dependent lower limit on the PQ scale 



/pq 



> 8.4 x 10 7 GeV 



VlGeW o*(T CT )i/4- ^ 



Almost the same limit applies to thermally produced sax- 
ions as well since their production is efficient only at high 
temperatures not far below Tr and leads basically to a 
thermal spectrum, i.e., (j3"3")l applies after substituting Td 
by T R and g* s (T D ) by 5 * S (T R ) ~ 228.75. 

Note that the saxions decay while being decoupled 
from the primordial plasma if T a <C Td or equivalcntly 



/pq 

-.1/3 



> 7.1 x 10 7 GeV i 



100 GeV/ 



^1/2 


' 232.5 " 




_g*{Ta)_ 



1/12 



(35) 



4 This g t g value accounts for the MSSM and the axion multiplet 
fields, which can all be considered as relativistic at Td if not only 
m a but also the axino mass satisfies raj <C Trj . 



If this condition is satisfied the axions emitted in those 
decays will not be thermalized and be free-streaming. 
Thus, the temperature T nr at which the non-thermally 
produced axions become non-relativistic reads 



T nr = 0.15 a; eV 



3.91 



10 TeV 



1/2 /io 9 GeV 



/pq 



g*s(T n 



1/3 



g*s{T a 



,1/3 
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(36) 



when defined via p a (T nr ) = m a . This shows that the 
emitted axions are expected to be still relativistic at the 
last scattering surface and even well thereafter for m a < 
10 TeV and x = 0(1). Thereby they can contribute to 
AA^ff even at late times where studies of the CMB and 
LSS allow us to probe the amount of radiation. 



A. BBN 

For Tj > 1 MeV, the axions from saxion decays con- 
tribute to the radiation density already at the onset of 
BBN and prior to e + e~ annihilation. This leads to a 
speed-up of the Hubble expansion rate and thereby to an 
output of 4 He that is more efficient than in standard BBN 
with AA^cff = 0. In turn, the inferred primordial 4 He 
abundance imposes upper limits on AN c s, whereas the 
inferred primordial D abundance constrains the baryon 
density Wb = Qhh 2 , with the normalized Hubble constant 
h ~ 0.7. Notably, two recent studies of the primordial 
4 He mass fraction Y p even report values in excess by more 
than 2cr over the standard BBN prediction: Izotov and 
Thuan HI find = 0.2565±0.001(stat.)±0.005(syst.) 
and Aver et at [47] V p Av = 0.2561±0.0108, with all errors 
refering to 68% intervals. As mentioned in the Introduc- 
tion, these results may be hints towards extra radiation 
at the onset of BBN, which can reside in the form of 
axions from decays of saxions from thermal processes. 

Based on recent studies of the primordial 4 He and D 
abundances [ill |4?], [66| and the recent Particle Data 
Group (PDG) recommendation for the free neutron life- 
time, r n = 880.1 ± 1.1 s [60| . we now derive AN e g limits 
from a BBN likelihood analysis, similar to the one in 
Ref. [67J , and explore the implications for the considered 
SUSY axion models. Relying on the Y p results given 
above and on the primordial D abundance reported by 
Petini et al. [H, log[D/H] p = -4.56 ± 0.04, we consider 
the two log-likelihood functions 



lnZ 



4 Ho 



1 (y p -0.2565r 

2 0.0051 2 
1 (y p - 0.2561) 2 

" 2 



lnT D = - 



for [46] 

0.0108^ for I 47 ! 

l(log[D/H] p T4.56) 2 
2 



(37) 
(38) 



0.04 2 

where small uncertainties related to nuclear reaction 
rates and also the ones related to the free neutron lifetime 
T n are not taken into account. Theoretical values for the 
primordial 4 He and D abundances are calculated with 



TABLE II. Constraints on AN c g from BBN and precision cos- 
mology. Based on the indicated data sets, the first two lines 
give the posterior maximum (p.m.) and the minimal 99.7% 
credible interval imposed by BBN using the prior AN c s > 
and after marginalization over uih- The third line quotes the 
mean and the 95% CL upper limit on AN e f! as obtained in the 
precision cosmology study of Ref. [4£| based on CMB data, the 
Sloan Digital Sky Survey (SDSS) data-release 7 halo power 
spectrum (HPS), and data from the Hubble Space Telescope 
(HST). 



Data 


p.m. /mean 


upper limit 


F P IT [46] + [D/H] p [66] 


0.76 


< 1.97 (3a) 


Y p Av [47] + [D/H] p [66] 


0.77 


< 3.53 (3<x) 


CMB + HPS + HST [48] 


1.73 


< 3.59 (2a) 



the BBN code PArthENoPE [68j using r„ = 880.1 s [60[ 
and < AN cff < 4 and 0.01 < u h < 0.03 as flat priors. 
Calculating the respective combined likelihood and after 
marginalizing over Wb, we obtain for AN c g the maximum 
likelihood posteriors and the minimal 99.7% credible in- 
tervals listed in the first two lines of Table ITTPI 

Let us now apply these BBN constraints to the case 
of extra radiation from saxion decays into axions. We 
evaluate AN eS (T) from ((3TJ) and ([321) for x = 1 and at 
T ~ 1 MeV, i.e., at the onset of BBN and above the tem- 
perature at which neutrinos decouple. Figure |4] shows the 
resulting AN e g contours in the m^-fpQ parameter plane 
as black (gray) lines for T R = 10 s (10 10 ) GeV. The solid 
(dashed) curves show - as labeled - the posterior max- 
imum (A7V ff)^'™j T ^ = 0.77 (0.76) and the upper limit 

(AiV eff )^ (IT) = 3.53(1.97), which disfavors the consid- 
ered region to its left by more than 3tr. The dotted lines 
indicate T„ — 1 and 10 MeV. The parameter region with 
T CT < 1 MeV is not considered since our BBN constraints 
on AiVoff do not apply to later decays@ Moreover, as de- 
scribed in the Introduction, additional cosmological con- 
straints can occur for T a < 1 MeV, which are beyond the 
scope of this work. 

For T a > 1 MeV, one sees that the BBN constraints 
on AiV e ff can disfavor significant regions of the m a -/pQ 
parameter plane in high Tr scenarios. These regions will 
become larger and move towards larger m a if x is smaller 



5 Note that the PDG-recommended value for the free neutron 
lifetime has changed recently from r n = 885.7 ± 0.8 s [69t to 
r n = 880.1 ± 1.1 s [§(j. If we use r n = 885.7 s in PArthENoPE, we 
can reproduce the posterior maxima and the minimal 95% cred- 
ible intervals given in the first two lines of Table III in Ref. [70l . 
In comparison, those posterior maxima are about 10% below the 
values obtained with T n = 880.1 s given in our Table HT1 

6 The calculations of PArthENoPE start at T = 10 MeV with 
the given AN c g values as input already at that temperature. 
The A7V c ff limits derived above are thus strictly applicable for 
T a > 10 MeV only. Modifications in the PArthENoPE code 
that will allow us to describe more accurately the region with 
T a < 10 MeV are postponed to future work. 




m CT [GeV] 



FIG. 4. Contours of AV cff at T ~ 1 MeV provided by non- 
thermally produced axions from decays of thermal saxions for 
x = 1, T ff > 1 MeV, and T R = 10 s GeV (black) and 10 10 GeV 
(gray). The BBN results given in Table ITT1 are illustrated by 
the solid (dashed) curves, which indicate the posterior max- 
imum (AV c ff)^'™j T j = 0.77 (0.76) and the upper 3a limit 
(AV e ff)^^ IT j = 3.53 (1.97) that disfavors the region to its 
left. The dotted lines show T CT = 1 and 10 MeV. 



than one but still sufficiently sizable such that the sax- 
ion decay into axions remains to be the dominant decay 
channel that governs r a . Moreover, the shown posterior- 
maxima contours illustrate that non-thermally produced 
axions from decays of thermal saxions can explain the ex- 
istence of extra radiation, AiV e g ~ 1, in agreement with 
the hints from BBN studies. For Tr > To, the shape of 
the AN e g contours is described by (|3"Tj) from decays of 
thermal relic saxions. The kink of the AN e g contours in- 
dicates the respective /pq value at which Xr = To- For 
larger /pq, Tr < Tb and (|32|) applies which is reflected 
by the Tr dependence of AN e g provided by axions from 
decays of thermally produced saxions. 

As mentioned in the Introduction, the energy density 
of the early Universe can receive contributions not only 
from thermal saxions but also from coherent oscillations 
of the saxion field [8j4l2|. In fact, axions from the de- 
cay of those non-thermal saxions constitute additional 
radiation as well [ll|> [l2j and can thereby increase AN e g 
already at the onset of BBN. However, in the parameter 
region considered in Fig. HI their contribution to AiVeg 
is basically negligible for an initial displacement of the 
saxion field from the vacuum of Oi ~ /pq- This can be 
seen, for example, in Fig. 1 of Ref. (TlJ where the energy 
density of coherent saxion oscillations is compared with 
an estimate of the one of thermal saxions. 

At this point, it is important to stress that significant 
additional restrictions are possible that depend on the 
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mass spectrum and on other aspects of the specific SUSY 
model possibly realized in nature. Here the masses of the 
axino and the gravitino are of particular importance since 
their thermal production can be very efficient in high Tr 
scenarios such as those explored in Fig. [4] 

For example, in R-parity conserving scenarios in which 
the gravitino is the lightest SUSY particle (LSP), the 
thermally produced gravitino density is limited from 
above by the dark matter density parameter wcdm = 
n CBM h 2 ~ 0.1 [60]. This disfavors T R > 10 10 GeV and 
translates into m 3/2 < 10 GeV for T R = 10 s GeV and 
universal gaugino masses at the scale of grand unification 
of m 1/2 - 0.5 TeV; cf. Fig. 2 in Ref. j^. For m a ~ m 3 / 2 , 
as expected in gravity-mediated SUSY breaking, this cos- 
mological constraint will then challenge the AJV e g ~ 1 
explanation presented in Fig. |4] for Tr = 10 s GeV and 
disfavor the one for Tr = 10 10 GeV. Depending on the 
next-to- lightest SUSY particle (NLSP), even more re- 
strictive upper limits on Tr are possible; cf. [37j and ref- 
erences therein. Additional limits related to axino cos- 
mology can be evaded, e.g., for ma > 2 TeV, a gluino 
mass of m s ~ 1 TeV, and / PQ - 10 10 GeV. Although 
possibly somewhat contrived from the model building 
point of view, the heavy axinos then decay typically 
before the NLSP freeze-out and the emitted sparticles 
will be thermalized such that the constraints associated 
with the NLSP will not be tightened. If a sneutrino is 
the NLSP (for which the NLSP-related constraints are 
rather mild (7ll-l73|). the shown AN e g ~ 1 explanation 
for Tr = 10 8 GeV can thereby turn out to be viable for 
Too- ~ ~ 10 GeV, where cold dark matter can reside 
in thermally produced gravitinos. 

In the alternative axino LSP case, one often finds more 
restrictive Tr constraints imposed by the dark matter 
constraint [23M25I l74j and also additional fpo constraints 
depending on the properties of the NLSP [74l476l | . Inter- 
estingly, these Tr constraints can be avoided in the case 
of alight axino LSP with m 5 < 0.2 keV (cf. Fig. 6 in [Hj]). 
Moreover, additional Tr constraints from BBN- imposed 
limits on hadronic/electromagnetic energy injection from 
late decayi ng g ravitinos can be evaded if the gravitino is 
the NLSP jH, [73. In such a setting, the lifetime of 
the gravitino NLSP is r 3/2 - 10 9 s (10 2 GeV/m 3/2 ) 3 and 
governed by its decay into the axino LSP and an ax- 
ion [22], [2} • While both of which are too weakly interact- 
ing to reprocess primordial nuclei, the emitted particles 
can contribute to AN c g at cosmic times t > r 3 / 2 [HI- 
Upper limits on AN e g imposed by CMB + LSS stud- 
ies have thereby been found to imply Tr < 10 11 GeV at 
the 5cr level for m 3/2 = 100 GeV and m s ~ 1 TeV [z|. 
This limit can be overly conservative since it does not 
include AN c g from saxion decays into axions. For rrig ~ 
1 TeV, which is still allowed by the ongoing LHC spar- 
ticle searches, m 3/2 - 100 GeV, and Tr - 10 10 GeV, 
gravitino decays into axions and axinos have been found 
to lead to A-/Vpff ~ 0.6 but only at times well after the 
BBN epoch [6l|, uM, ■ Taking into account the additional 
contribution to AN e f[ at such late times from saxion de- 



cays (which we consider explicitly below), we find that 
this point in parameter space remains to be allowed. 
This implies viability of the corresponding explanation 
of AN cS ~ 1 at the onset of BBN for Tr = 10 10 GeV 
shown in Fig. |4j Here one can easily accommodate also 
the small additional contribution of AN c g < 0.017 pro- 
vided by light thermal axinos at the onset of BBN |76| • 
For fpQ ~ 10 12 GeV, cold dark matter can then reside 
in the form of an axion condensate (cf. Fig. [7] below) 
whereas axinos will be hot dark matter [24[ with associ- 
ated LSS constraints imposing mj < 37 eV [76| . While 
the lightest ordinary sparticle (LOSP) can still be long 
lived, BBN constraints related to its decay can be evaded. 
For the stau LOSP case, this is illustrated explicitly in 
Fig. 21 of Ref. [76j|. Thereby, one arrives at viable sce- 
narios with different AN e g predictions at the onset of 
BBN and much later, in which even Tr ~ 10 10 GeV is 
possible, e.g., allowing for the explanation of the baryon 
asymmetry via thermal leptogenesis (7{| . 

B. CMB and LSS 

Axions from saxion decays can contribute to AN e g at 
the CMB decoupling epoch, even for T CT > 1 MeV, as 
described below (f3"6"| . Extra radiation at that epoch de- 
lays the time of radiation-matter equality and is probed 
by studies of the CMB anisotropics and the LSS distri- 
bution. Here hints towards N c s > 3 have been found 
that are more pronounced than those from BBN con- 
sidered above; see [48j - [52| and references therein. For 
example, the Wilkinson Microwave Anisotropy Probe 
(WMAP) collaboration finds a 68% credible interval of 
N c ff = 4.34+Qjg [5l[ when combining their 7-year data 
with measurements of the baryonic acoustic oscillation 
(BAO) scale and todays Hubble constant Hq. Another 
precision cosmology stud y a rrives at a 95% credible in- 
terval of N cS = 4.781^75 when combining CMB data 
with data from the Sloan Digital Sky Survey data-release 
7 halo power spectrum (HPS) and the Hubble Space Tele- 
scope (HST). Based on this combined CMB + HPS + 
HST data set, we use the mean for AN c g and the 95% 
CL upper limit on AN c g, as quoted in TableQll to explore 
implications for the considered SUSY axion models. 

Evaluating AN ca (T) from $tt$ and (j3"2"|) for x = 1 
and at T <C 1 MeVQ we obtain the AN c g contours 
for T R = 10 8 (10 10 ) GeV as shown by the black (gray) 
lines in Fig. [5j The solid lines indicate the upper limit 
AiVeff = 3.59 with the to ct _ /pq parameter regions to their 
left disfavored at the 2a level by the CMB + HPS + HST 
data set. The dashed lines show the corresponding mean 



7 Note that our theoretical results for AA r of j(T) at T ~ 1 MeV and 
at T < 1 MeV agree. The T dependence in ((3TJ) and ((32j results 
from the factor (T/T„) 4 [g» s (T)/10.75] 4 / 3 . This factor equals 
one for T ~ 1 MeV, where T v = T and g*s(T) = 10.75, and for 
T < 1 MeV, where T„ = (11/4)" 1 / 3 T and g» S \T) = 3.91. 
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FIG. 5. Contours of AiV cff at T < 1 MeV provided by non- 
thermally produced axions from decays of thermal saxions 
for x = 1 and Tr = 10 s GeV (black) and 10 10 GeV (gray). 
The solid curve shows the 2a limit AN e a = 3.59 and the 
dashed curve the mean AN e « = 1.73 based on the CMB + 
HPS + HST result [U quoted in Table |TI] The dash-dotted 
curve indicates AN c g = 0.26 which is the expected 68% CL 
sensitivity of the Planck satellite mission p3l , |54| . On the 
dotted lines, T a = l and 10 MeV, as in Fig. H 
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FIG. 6. Contours of AiV off at T < 1 MeV provided by non- 
thermally produced axions from decays of thermal saxions for 
x = 1 and f PQ = 10 10 GeV (light gray), 10 11 GeV (black), and 
10 12 GeV (gray). The solid lines show the 2a CMB + HPS + 
HST constraint AN c g < 3.59 which imposes upper limits on 
the reheating temperature Tr. For /pq = 10 10 GeV, the Tr 
limit appears at smaller m a outside of the considered range. 
The dash-dotted contours indicate the expected Planck sen- 
sitivity of AiVeff = 0.26. 



AiV eff = 1.73 and the dash-dotted lines AAT off = 0.26. 
The latter is the expected 68% CL accuracy of the Planck 
satellite mission [53, [54j mentioned already in the Intro- 
duction. To guide the eye, we show again the T a — 1 and 
10 MeV contours as dotted lines. Here we can provide 
the AN e g contours also in the region with T a — 1 MeV. 
However, as described in the Introduction, additional re- 
strictive constraints are expected in that region. 

In comparison to the BBN-imposed limits, one finds 
that the shown 2a limit from precision cosmology disfa- 
vors basically the same parameter region as the conser- 
vative 3cr limit shown in Fig. 01 The results from preci- 
sion cosmology thereby allow for a slightly larger AN e g 
at given values of m CT , /pq, and Tr. Indeed, comparing 
the mean and the 2a limit in Table lU with the posterior 
maxima and the 3er limits from the BBN study, one finds 
a potential hint towards a AN c g value at T -C 1 MeV 
that is higher than the one at T ~ 1 MeV. As discussed 
already in the preceding section, this may be a first hint 
for AN e g ~ 1 already prior to BBN due to axions from 
saxion decays plus an additional late contribution from 
gravitino NLSP decays into axions and LSP axinos such 
that AN e g ~ 2 at the time of CMB recombination. With 
the expected A N e $ sensitivity of the Planck satellite mis- 
sion, this possibility will be tested further soon. More- 
over, for scenarios in which axions from saxion decays are 
the only significant source for AN e s, the Planck results 
will allow us to probe significant regions of the m l7 -/pQ 



parameter space which have not been accessible by AN c h 
studies so far. For fixed Tr values, this is indicated by 
the dot-dashed lines in Fig. [5] 

The limits shown in Figs. |4] and [5] in the t?t,<t _ /pq plane 
for fixed Tr values can be translated into upper limits on 
the reheating temperature Ta. In Fig. [5] the solid lines 
show the upper limits on Tr imposed by the 2a CMB 
+ HPS + HST constraint AAc ff < 3.59 as a function of 
TTtfj for x = 1 and / PQ = 10 11 GeV (black) and 10 12 GeV 
(gray). The expected Planck sensitivity AN e g = 0.26 
is indicated by the corresponding dash-dotted lines and 
in light gray for /pq = 10 10 GeV. Note that the upper 
limit does not show up for the latter /pq value in the 
considered m a range. The Tr dependence of the con- 
tours is described by (j32|) and disappears for cosmolog- 
ical scenarios with Tr ^ Td where (|3"Tj) applies. Here 
we should stress that the shown upper limits on Tr rely 
on non-thermally produced axions from decays of ther- 
mal saxions providing the only significant contribution to 
AA^ff at T <C 1 MeV. In scenarios with additional siz- 
able contributions (e.g., from late gravitino NLSP decays 
into an axino LSP and the axion), the considered extra 
radiation constraint will impose more restrictive Tr lim- 
its. Nevertheless, the shown upper limits will remain to 
be applicable as conservative guaranteed limits. 

Let us compare our results shown in Fig. [S] with exist- 
ing results. For example, the yellow curve in Fig. 5(a) 
of Ref. presents an upper Tr limit imposed by 
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AN eS < 1 that disfavors basically T R > 10 6 GeV for 
/pq = 10 10 GeV and the whole m a range considered 
above. With the assumed initial saxion field displace- 
ment of Ci ~ /pq, that existing limit is governed also by 
thermal saxions that decay into axions. However, we find 
that it is overly restrictive due to the omission of the fac- 
tor [g*s{Ta)/g* s {T)} 4 / 3 in Eq. (24) of Ref. 0. Thereby 
our AJVeff expression (|30p shows different dependences on 
g*s{T) and g*s(T a ). Remaining differences are due to the 
result from our explicit calculation of the thermal saxion 
production rate and the different definitions of the PQ 
scale addressed already in footnotes 2 and 3 above. As 
a result, we find that a large part of the (m a , /pq, Tr) 
region previously thought to be excluded is actually not 
restricted by the amount of additional radiation from late 
decays of thermal saxions. 



VIII. RELIC AXION DENSITY 

We find it instructive to compare the density param- 
eters of three different axion populations that can be 
present today in the considered SUSY axion models: 
(i) Qa q ^ TF h 2 of thermal relic/thermally produced ax- 
ions, (ii) fl^ TP h 2 of non-thermally produced axions from 
decays of thermal saxions, and (iii) f2„ is h 2 of the ax- 
ion condensate from the misalignment mechanism. The 
latter originates from coherent oscillations of the axion 
field after it acquires a mass due to instanton effects at 
T < 1 GeV. This is the axion population that can pro- 
vide the cold dark matter in our Universe, as mentioned 
at the end of Sect. lVILAl For details on this misalignment 
mechanism we refer to 0, [H, 0] and references therein. 
Here we quote the density parameter, 



0.15 0? 



10 12 GeV 



(39) 



which is governed by the initial misalignment angle 9i 
of the axion field. This expression applies to non-SUSY 
and SUSY settings. In the considered case in which the 
PQ symmetry breaks before inflation and is not restored 
thereafter, Tr < /pq, a single 6i value will enter (|39| . 
Since thermal processes involving axions proceed at neg- 
ligible rates at T < 1 GeV for /pq respecting ([5]), the 
axion condensate cannot be thermalized. Accordingly, 
nf is h 2 can coexist with n c a q/TP h 2 and ^ TP /i 2 , which 
we calculate in the following. 

Since thermal relic and thermally produced axions 
have (basically) a thermal spectrum, one can describe 
the associated density parameter approximately by 



n e q/T P^2 ^ jipf o) 2 +m 2 a ^cq/TP s{T ^ , p ^ (40 ) 

where p c /[s(T )h 2 } = 3.6 eV and V a cq/TP = V CT cq/TP as 
described in Sect. |V] With the present CMB tempera- 
ture of T = 0.235 meV and an axion temperature to- 
day of T a , = [.g* s (7o)/228.75] 1 / 3 r ~ 0.06 meV, the 



average momentum of thermal axions today is given by 
(Poo) = 2.701 T a fi. When comparing this momentum 
with the axion mass m a , one finds that this axion pop- 
ulation is still relativistic today for /pq > 10 11 GeV. At 
and before the CMB decoupling epoch, T > 1 eV, axions 
from thermal processes were relativistic for /pq in the 
full allowed range ©. In the considered SUSY settings, 
they contribute at most 



AiV off (T) 
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'g*s{T)~ 


4/3 
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_228.75_ 







, (41) 



i.e., AN cS (T) < 0.0097 for 10 MeV > T > 1 eV, which is 
far below the Planck sensitivity and easily accommodated 
by the AN e g limits discussed above. 

The density parameter of non-thermal axions emitted 
in late decays of saxions from thermal processes reads 



n 



NTP.2 



+ mi Y. 



2 V eq/TP 



s(T )h 2 /p c (42) 



with the present momentum of these axions given by 



NTP 
Pa,0 



g*s{Tg) 
g*s{T a ) 



"I 1/3 



To 



(43) 



when applying the sudden decay approximation. Thus, 
n^ TP h 2 will depend on m a if these axions are still rel- 
ativistic today, i.e., when To > T nr given by (|3"6")l above. 
As extensively discussed in the previous section, this non- 
thermal axion population can provide a significant con- 
tribution to AN e ff prior to BBN and thereafter. In fact, 
to relate n^ TP h 2 to this AiV off : 



one can use 



n 



NTP.2 



h z = { 1.1 x 10" 



10 9 GeV \ 2 
/pq ) { ID" 3 



+ 3.2 x 10 



-ii 



10.75 



g*s(T) 



8/3 ~) 

A7V e 2 ff (T) , (44) 



where the m a dependence is now absorbed into AN e g(T). 
Thus, the discussed AiV ff constraints translate directly 
into upper limits on fl^ TP h 2 . For T < 1 MeV and /pq 
such that the first term on the right-hand side of (|4"4")) 
is negligible, those limits are described by fl^ TP h 2 = 
5.7 x 10~ 6 AN e g. (This applies to thermal axions as well 
if they are still relativistic today.) 

Figure [7] shows Q^ TP h 2 contours that correspond to 
AV cff values of 3.59 (solid), 1.73 (dashed), and 0.26 
(dash-dotted). As in Fig. [5j these values are motivated 
by the CMB + HPS + HST result 0] quoted in TableU 
and the expected 68% CL sensitivity of the Planck satel- 
lite mission [53|, Contours of Qa q ^ TP h 2 are shown 
for T R = 10 8 (10 10 ) GeV by the solid black (gray) lines 
and for larger Tr > Td by the unlabeled dotted line. 
On this dotted line at / PQ > 10 11 GeV, (gT|) applies so 
that AiV cfT = 0.0097 and Sl^h 2 = 5.5 x 10" 8 reside in 
thermal relic axions that are still relativistic today. The 
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Taking into account the relation between /pq and 
vn a , the analog of a Lee- Weinberg curve is given by 



10 10 
/pq [GeV] 

FIG. 7. The density parameters of the axion condensate 
from the misalignment mechanism Q^ is h 2 for 8i = 0.01, 
0.1, and 1 (dotted lines), of non-thermally produced axions 
from decays of thermal saxions Q% TF h 2 for AN cS = 3.59 
(solid), 1.73 (dashed), and 0.26 (dash-dotted), and of ther- 
mal relic /thermally produced axions Q'^ 1 ^ TP h 2 for Tr = 10 8 
(black) and 10 10 GeV (gray). The dotted line connected to 
the latter indicates Q, c ^h 2 for larger Tr with Tr > T D . The 
dark matter density parameter Clcnmh 2 ~ 0.1 60J is indi- 
cated by the horizontal gray bar. As in Fig. [5] we consider 
AN cB at T < 1 and show values based on the CMB + HPS + 
HST result [S| quoted in Table [TT] and the expected 68% CL 
sensitivity of the Planck satellite mission [Hj • 



n 



NTP u2 



^eq/TP^ 2 anc j can j n _ 



ferred from Fig. [7J Depending on the initial displacement 
of the saxion field from the vacuum, <Ji, and on the mass 
spectrum, there can be additional contributions to the 
axion density parameter, e.g., from decays of the saxion 
condensate into axions and/or a gravitino NLSP into ax- 
ions and LSP axinos. In such cases, sizable additional 
contributions also to AN c g are possible which will affect 
the Sl^ TP /i 2 contours in Fig. [7J Thus the shown contours 
should be understood as conservative maximum values. 

One can consider Fig. [7] as a SUSY generalization of 
Fig. 4 in Ref. 17], which allows one to infer the axion 
analog of the Lee- Weinberg curve in non-SUSY scenar- 
ios. Whereas fl^ TP h 2 can govern the axion density for 
small 9i and/or small /pq in non-SUSY scenarios [l7l |. 

we find nfj p h 2 > 2fC /TF V in the considered SUSY 
scenarios. This can be seen in Fig. [Jj and when com- 
paring (|4TJ|) and (j4"2"]l . If SUSY and a hadronic axion 
model are realized in nature, the axion density param- 
eter can thus be governed by non-thermal axions from 
decays of thermal saxions and/or the axion condensate 
from the misalignment mechanism. Interestingly, both 
of these populations may be accessible experimentally: 
While signals of the axion condensate are expected in 
direct axion dark matter searches [8(|, the findings of 
AiV e )f studies may already be first hints for the existence 
of non-thermal axions from saxion decays. 



IX. CONCLUSION 



labeled dotted lines indicate Sl^ IS /i 2 of the axion con- 
densate from the misalignment mechanism for 6i = 0.01, 
0.1, and 1. For 0* ~ 1 and / PQ - 10 12 GeV, this cold 
axion population can explain the dark matter density 
Ocdm^ 2 — 0.1 [60] displayed by the gray bar. 

Considering the 2a limit AY c ff < 3.59 in Fig. [7J one 
sees that it constrains f2^ TP /i 2 to values that stay below 
the photon density fi 7 /i 2 ~ 2.5 x 10~ 5 [601 ] . Remarkably, 
Planck results are expected to probe even much smaller 
f^ TP . The testable values can be as small as an or- 
der of magnitude below f2 7 if axions emitted in decays 
of thermal saxions are the only significant contribution 
to AiVpff. In contrast and similarly to the non-SUSY 
case [13|, it will remain to be extremely challenging to 
probe the axion population from thermal processes with 
its small contribution of AN e g < 0.01. 

Note that m a changes along the n^ TP h 2 curves in 
Fig.UJfor fixed Tr, and x since we indicate results for fixed 
values of AY c ff . Indeed, additional BBN constraints can 
disfavor parts of the shown contours when T a < 1 MeV. 
For T a > 1 MeV, BBN constraints on AjV eff - such as 
the ones considered in Fig. 0] - can also be displayed in 
terms of il^ TP h 2 . On the logarithmic scale considered in 
Fig. [JJ they are similar to the shown ones. 



We have explored thermal production processes of ax- 
ions and saxions in the primordial plasma, resulting ax- 
ion populations and their manifestations in the form of 
extra radiation AiV e ff prior to BBN and well thereafter. 
The considered SUSY axion models are attractive for a 
number of reasons. For example, they allow for simulta- 
neous solutions of the strong CP problem, the hierarchy 
problem, and the dark matter problem. 

Here we have focussed on the saxion, which can be a 
late decaying particle and as such be subject to various 
cosmological constraints. We find that the saxion decay 
into two axions is often the dominating one. For a saxion 
mass of m a > 1 GeV, such decays occur typically before 
the onset of BBN. We have shown that the emitted ax- 
ions can then still be relativistic at the large scattering 
surface. Thereby, they can provide sizable contributions 
to extra radiation A7V c ff that is testable in BBN studies 
and in precision cosmology of the CMB and the LSS. 

We have aimed at a consistent description of both the 
thermal axion/saxion production and of saxion decays 
into axions. This has motivated our careful derivations of 
the Lagrangian £pg that describes the interactions of the 
PQ multiplet with quarks, gluons, squarks, and gluinos 
and of £p,2, that describes the interactions of saxions 
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with axions in addition to their kinetic terms. The re- 
quirement of canonically normalized kinetic terms defines 
the scale vpq , which governs the saxion-axion-interaction 
strength together with another PQ-model-dependent pa- 
rameter x < 1. On the other hand, the form of the ef- 
fective axion-gluon-interaction term defines the PQ scale 
/pq. Considering the emergence of this term from loops 
of heavy KSVZ fields in an explicit hadronic axion model, 
we find fpQ = \/2vpq. This is in contrast to numerous 
existing studies which treat vpq and /pq synonymously. 

Relying on the derived form of CpQ , we have calculated 
the thermal production rates of saxions and axions and 
the resulting yields in the hot early Universe. Despite 
differences in the interaction terms, we find that the rate 
for thermal saxion production agrees with the one for 
thermal axion production. This implies an agreement 
also of the calculated thermally produced yields and of 
our estimates of the decoup ling temperatures Td- By 
applying HTL resummation 30J and the Braaten-Yuan 
prescription [3l|, finite results are obtained in a gauge- 
invariant way consistent to leading order in the coupling 
constant and screening effects are treated systematically. 

Using our result for the thermally produced saxion 
yield, we have calculated AN e s provided in the form of 
axions from decays of thermal saxions. This has allowed 
us to demonstrate that such a AN c g contribution can in- 
deed explain the trends towards extra radiation beyond 
the SM seen in recent studies of BBN, CMB, and LSS. 

To account for the current PDG recommendation for 
the free neutron lifetime, r n = 880.1 ±1.1 s [60l |. 
we have performed a BBN likelihood analysis with 
PArthENoPE [68[ and based on recent insi ghts on the pri- 
mordial abundances of 4 He [H|47| and D |66j]. For A7V eff 
at the onset of BBN, we thereby obtain posterior max- 
ima of 0.76 and 0.77 and 3<r upper limits of 1.97 and 3.53 
with the Y p results of [46| and 47 1, respectively. When 



comparing these values with results from studies of the 
CMB and LSS, we find that the latter provide compat- 
ible but more pronounced hints for extra radiation. For 
example, the precision cosmology study of [48[ reports a 
mean of 1.73 and a 2a limit 3.59 for AN cS at T < 1 MeV 
when using the CMB + HPS + HST data set. 

We have translated the upper limits on AN e g quoted 
above into bounds on /pq, m CT , and Tr. These bounds 
can disfavor significant regions of the m CT -/pQ parameter 
plane in high Tr scenarios. However, we find that our 
limits leave open a considerable parameter region previ- 
ously thought to be excluded [ll| . Significant parts of the 
allowed parameter region have been identified, which will 
become accessible very soon with the upcoming results 
from the Planck satellite mission. 

The explanation of the above hints for extra radiation 
via axions from decays of thermal saxions requires a rel- 
atively high reheating temperature of Tr > 10 7 GeV for 
ni a > 0.1 GeV. Such high Tr scenarios can be in conflict 
with cosmological constraints due to overly efficient ther- 
mal production of axinos and gravitinos. To illustrate the 



viability of AN e s ~ 1 from saxion decays, we have de- 
scribed two exemplary SUSY scenarios which allow for 
T R = 10 s and 10 10 GeV: 

(i) With the gravitino LSP as cold dark matter and a 
sneutrino NLSP, the presented AN c g ~ 1 explanation 
for Tr = 10 s GeV can be viable for m a ~ TO3/2 ~ 
10 GeV and rrig ~ 1 TeV. This explanation requires 
/pq ~ 10 10 GeV and heavy axinos, nia > 2 TeV, which 
decay prior to NLSP decoupling. Here AN c g ~ 1 is al- 
ready present at the onset of BBN and does not change 
thereafter. Accordingly, we expect that the Planck re- 
sults will point to a AN e g value that is consistent with 
the one inferred from BBN studies. 

(ii) With a very light axino LSP, ma < 37 eV, as hot 
dark matter and a gravitino NLSP, the AN c g ~ 1 ex- 
planation for Tr = 10 10 GeV can be viable for m a ~ 
m 3/2 ~ 100 GeV and mg ~ 1 TeV. Here this explana- 
tion requires /pq ~ 10 12 GeV so that cold dark matter 
can be provided by the axion misalignment mechanism. 
With the stau as the LOSP, further potential BBN con- 
straints can be evaded. Note that Tr = 10 10 GeV allows 
for successful thermal leptogenesis. The saxion decays 
give AiVcff ~ 1 already at the onset of BBN. However, 
late gravitino NLSP decays into the axion and the axino 
LSP can provide an additional contribution of AN c g ~ l 
well after BBN [H Izl- Thus, it will be interesting to 
see whether the Planck results confirm the trend towards 
an excess of extra radiation that is more pronounced at 
late times. For example, the finding of AN e g ~ 2 at late 
times will be a possible signature expected in this setting. 

If a SUSY hadronic axion model is realized in nature, 
three different axion populations will be present today: 
thermally produced/thermal relic axions, non-thermally 
produced axions from decays of thermal saxions, and the 
axion condensate from the misalignment mechanism. We 
have calculated and compared the associated density pa- 
rameters. The results allow us to infer the axion analog 
of the Lee-Weinberg curve. For /pq > 10 12 GeV and 
an initial misalignment angle of 9i ~ 1, the axion den- 
sity parameter is governed by the axion condensate. In 
that parameter region this population may be accessible 
in direct axion dark matter searches. For smaller /pq 
and smaller 0i, axions from saxion decays can dominate 
the axion density parameter. While it will be extremely 
challenging to probe thermal axions, Planck may confirm 
AN e g signals of this non-thermally produced population 
in the full allowed /pq range. Since the considered ax- 
ion populations can coexist, there is the exciting chance 
to see signals of both axion dark matter and axion dark 
radiation in current and future experiments. 
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